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ABSTRACT 


The cometary bow wave of P/Glacoblnl-Zlnner 
has been analyzed using the complete set of ICE 
field and particle observations to determine If 
It Is a shock. Changes In the magnetic field 
and plasma flow velocities from upstream to 
downstream have been analyzed to determine the 
direction of the normal and the propagation 
velocity of the bow wave. The velocity has 
then been compared with the fast magnetosonlc 
wave speed upstream to derive the Mach number 
and establish whether It Is "supersonic", i.e., 
a shock, or "subsonic," i.e., a large amplitude 
•wave. The various measurements have also been 
' compared with values derived from a Ranklne- 
Hugonlot analysis. The results Indicate that, 
Inbound, the bow wave is a shock with M ■ 1.5. 
Outbound, a subsonic mach number Is obtained, 
however, arguments are presented that the bow 
wave Is also likely to be a shock at this 
location. 


1. INTRODUCTION 

The analysis referred to In the title addresses 
cwo basic questions: (1) Is the bow wave a 
shock? (2) can the observed properties up-and 
downstream of the bow wave be reconciled with 
the Ranklne-Hugoniot relations? 

The basic approach to answering the first of 
these questions is to derive the velocity of 
the bow wave in the upstream solar wind and 
determine whether the Mach number is > 1 or < 
1, I.e., supersonic or subsonic. At this stage 
of the analysis, the reconciliation with the 
Ranklne-Hugoniot relations Is principally re- 
stricted to comparisons between observed and 
calculated jumps In the particle density 
and magnetic field strength. We use the mathe- 
matical formulation to be found In Tidman and 
Krall (Ref. 1, p. 11) which Is basically a sin- 
gle fluid description. 

2. ANALYSIS PROCEDURE 


SHOCK 

COORDINATES 



Figure 1. Shock coordinates. The X axis Is the 
shock normal and the Y axis Is In the 
direction of the Jump In the vector 
magnetic field. The other parameters, 
which lie in the X-Y plane, are defined 
In the text. 

normal to the shock surface. As Is customary, 
we make use of the coplanarlty of the velocity 
vectors, Vj , V 2 , trie magnetic field vectors Fj , 

8*2 • and the normal, ft. A representation of 
these vectors appears In Figure 1 In coordin- 
ates aligned with, and moving with, the shock. 

The principal axes are (1) the normal, ft, where 
the symbol, a , Indicates a unit vector, (2) 
the vector change in the field from upstream to 
downstream, A B - F 2 - B^, and (3) the ortho- 


The analysis procedure that has been followed is 
based upon experience obtained In studying In- 
terplanetary shocks and planetary bow shocks 
(Ref. 2, 3). The first step Is to determine the 


gonal direction parallel tc the shock, also the 
direction of the electric field. The analysis 
Is based on three dimensional velocity vectors 
obtained by the Energetic Particle Anisotropy 
Spectrometer (EPAS) (Ref. 4). Three dimensional 
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measurements of anisotropies in the pick ur 
Ions are used to find a reference system movtnj 
1 with velocity, V t which renders the distribu- 
tions isotropic (Ref, 5). 

U’e use the formulation first Introduced by 
Abraham-Schrauner (Re_f. 6) which Involves the 
cross products of A V and A B as shown In che 
figure. Our experience with other shocks In- 
dicates that this approach Is generally reli- 
able even In che case of perpendicular shocks, 
l.e., when the angle, 0 Bn , between n and the 
upstream magnetic field, B^, Is a 90*. The 
computation of 0 Bn is the second step In our 
analysis. 

The speed of the bow wave In the solar wind 
frame, v, , la calculated ualng two complemen- 
tary methods. The more conventional approach 
Is based on che conservation of mass and can 
be expressed as Vg ■ I H 2 / ^ ^ 2 " N^)) A V • n 
where Hi, Hj are the fluid densities up-and 
downstream. The other method, which we deve- 
loped as an alternative and which depends on 
che magnetic field rather chan Che particle 
density, Is a condition satisfied by all kinds 
of electromagnetic waves. Basically, che wave 
velocity equals the rado of the wave electric 
field, in chls Instance simply A V x IT,, t;o che 
transverse magne tic field change associated 
with the wave , | A B|. An alternative deriva- 
tion of this relation Is based on recognizing 
that, In the shock frame, che wave Is non-pro- 
pagating and Che eleccrtc field is conserved 
( A t ■ 0) across che shock. The shock speed 
In Inertial space, Vj , Is also calculated by 
1 finding che component of che upstream solar 
''wind velocity along ft, l.e., Vj ■ Vg + Vi • ft. 

Before obtaining the Mach number, It Is neces- 
sary to derive che phase velocity for che fast 
mode of propagation upstream of che bow wave. 
The conventional equations were used for the 
Alfven speed, ■ B/(4 10 )*'^, che Ion sound 
speed, C s - ( y p/ p )*' , and the fast mode 
speed, Cf, which depends on C^, C s and che 
angle of propagation with respecc to che mag- 
nedc field, 6 Bn (Ref. 1, p. 15). The pres- 
sure, p, Is che sum of the partial pressures 
associated wlch the solar wind electrons, the 
solar wind protons and the heavy cometary Ions. 
The electron and heavy Ion properties were 
measured or Inferred from che ICE measurements 
(Ref. 7,8). Since no proton measurements are 
available, the wave speeds depend on reasonable 
assumptions for the proton densities and tem- 
perature. Ue have followed the policy of 
varying all the densities and cemperatures over 
reasonable limits In order to carry ouc a 
formal error analysis and ensure chat uncer- 
tainties In these quantities would not change 
our conclusions. The density, p , Is the sum of 
the density-mass products for the same three 
constituents. In the expression for Cj, y Is, 
of course, the ratio of specific heats for 
which we assume the conventional value of 5/3. 

Finally, the Mach number is calculated from H » 
V s /C f . This quantity essentially provides our 
definition of a shock, l.e., It Is a large 
amplitude wave propagating at supersonic 
speeds . 


UPSTREAM DOWNSTREAM 

|-H Hi 




z 

CO 

Ui 

oc 

3 

CO 

CO 

UJ 

c* 

&. 

mrJ 

< 


01 ’CO 


09; 15 


09i30 


09.45 


10O0 


Figure 2. Observations of the Inbound bow wave. 

Three of che principal pressures (mag- 
netic, solar wind electron and aolar 
wind convective pressures) are shown, 
The horizontal bars above the daca 
show the analysis Intervals correspond- 
ing to the results In Table I. 


An example of che daca used In the analysis Is 
shown In Figure 2. The quantities plotted are 
actually three of the principal pressures: the 
magnetic pressure, B^/8 x , the electron pres- 
sure, n e kT e , and the convective pressure of the 
solar wind, n e M p H, where M p Is che proton 
mass. The changes in these quantities across 
the Inbound bow wave are shown with the analy- 
sis Intervals used to obtain the results pre- 
sented below Indicated at che cop of the fi- 
gure. Overall, the decrease in convective 
pressure and the accompanying Increases In 
magnetic and internal pressure are evident In 
spice of the intense hydromagnetlc turbulence 
(Ref. 9) that Is present. 

The large irregularities Introduced by the tur- 
bulence pose a significant problem by making It 
difficult to determine accurately che para- 
meters needed In the analysis. In fact, we 
have found chat a variety of solutions can be 
generated by varying the location and/or dura- 
tion of the analysis Intervals. Other analysis 
procedures can also lead to alternative solu- 
tions for che same analysis Intervals, e.g., we 
have compared the analysis described above wlch 
the output of the program recently developed by 
Vlnas and Scudder (Ref. 10). It Is, therefore, 
Important to establish criteria to be met by an 
acceptable solution (or solutions). These cri- 
teria will be described below as we consider 
the results obtained so far. 



TABLE ! . 
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• n: 0.1210, -0.7089, -0.6948 


* e Bn : 8(U 


• V. 


AVXB, 


PERPENDICULAR 
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189 km /s 

• V s • (Ng/Ng-NjlAV-fi ■ 189 km/ s 

• Vj • n * -102, V| ■ 86 km/ s 

• C A - 68, C s - 106 km/s 




125 km/ : 


189 


• M = - 1.5 > 1.0 

125 

• A SHOCK 


3. ANALYSIS RESULTS 

Table I contains the parameters appropriate to 
the inbound crossing of the bow wave, The first 
row contains the components of ft in solar-eclip- 
tic coordinates (X toward the sun, Y parallel 
to the ecliptic and positive to the East, and Z 
parallel to the north ecliptic pole). One of 
our acceptability criteria Is that the normal 
should have a physically reasonable orienta- 
tion. Figure 3 shows the encounter geometry as 
seen looking toward the sun along the tall of 
Clacoblni-Zlnner. It can be seen that the 
Intercept Was basically from south to north so 
that a reasonable shock normal should have ny < 
0, n z < Cl Inbound and n y > 0, n z > 0 outbound. 
In both cases, n should have a component point- 
ing Into the upstream solar wind, so that n x > 
0 . 

The angle, 8 , is * 80* so that the wave or 
shock is nearly perpendicular as anticipated 
from the geometry of the trajectory and the 
orientation of the interplanetary magnetic 
field at the time of the crossing. 



Figure 3. Encounter Trajectory Looking Sun- 
ward. Contours corresponding to 
a bow shock and plasma tall are 
shown schematically along with the 
coordinate systen used in the 
analysis (upper right). The normals 
for the inbound and outbound cros- 
sings of the bow wave or shock are 
indicated . 


The values of and C s show that the Ion sound 
speed Is dominant. Their ratio implies an up- 
stream 8 - 8 a p/B 2 of 2.9. This value, which 
Is larger than the value near 1.0 that la typi- 
cally observed In the solar wind, demonstrates 
the Influence of the cometary pick-up Ions 
which have come to dominate the energy density 
of the plasma upstream of the bow wave. For 
this quasi-perpendicular orientation, the frst 
mode speed is essentially Cj « (C^ + C^ 2 ) 1 ' , 
leading to a value of 123 km/s. 


The two methods of computing V g lead to identi- 
cal results of 189 km/s. Although this exa<Yt 
agreement Is undoubtedly a coincidence In this 
case, we generally require acceptable solutions 
to yield approximately equal values of V 3 . The 
Inertial speed, * 86 km/s, Is not zero as 
would be expected for a bow wave or standing 
shock in an ideal situation. From the analysis 
of planetary bow shocks, however, it has been 
found that the Inertial speed 13 rarely zero 
but that the shocks are non-stationary with 
Instantaneous speeds of tens of km/s up to 
values of * 100 km/s. This realization forms 
the basis of another of our criteria, t.e., Vj 
need not be zero but should be small (much less 
than the solar wind speed). 


Hence, the Mach number, M ■ 189/125 ■ 1.5 > 1.0. 
The conclusion to be drawn Is that the Inbound 
crossing Is, In fact, a shock. Further support 
for this conclusion Is provided by the Ranklne- 
Hugonlot solutions, (with 8 ■ 2.9 and 8 - 80*) 
which show that the ratios N 2 /N 1 , ®2 / B 1 as 
observed are consistent with M ■ 1.6 + 0.1. 

The outbound crossing is considered next, the 
results showing that it is not easy to obtain 
unambiguous results In this case. Table 2 
contains the results for an acceptable solution. 
As above, we have insisted that ft have a rea- 
sonable orientation (n x >0, ny > 0, n z > 0 and 
demonstrating a reasonable quantitative rela- 



TABLE 1! . 


A PROBLEM OUTBOUND: 


• fij 0.4236, 0.1822, 0.8873 


FOR PERPENDICULAR SHOCK OR WAVE 


“W 87 - 7 


PERPENDICULAR 


v $ • |avxb 2 |ab - 60 km/ s | 


m ±13 


V s - (N 2 /N 2 -N 1 )AV*n«86 km/ sj 
Vj • n = 0, V ( ■ 60-86 km/ s 
C A ^ 32 km/ s, C s - 94 km/ s 


Cj - 99 km/ s 

• M f s |f - 0.73 +0.13 <1.0 

• SUBSONIC - A WAVE? 


tlonshlp), The wave turns out again to be near- 
ly perpendicular with 0 Bn " ®8*. tuo 
values of the shock speed are slightly different 
with V g » 73 + 13 km/s, the average value being 
low compared to that obtained Inbound. Again, 
the lnerclal speed of 60-86 km/s appears within 
an acceptable range. In the solar wind up- 
stream of the outbound crossing, , has de- 
creased to 32 km/s principally as a result of a 
decrease 'in the field magnitude by a factor of 
approximately 2. The sound speed Is also re- 
duced slightly so that Cf - 99 km/s. However, 
even this reduction Is not adequate to yield a 
Mach number greater than one. Thus, M “ 73/99 
- 0.73 + 0.13 < 1.0. 

Taken at face value, this result Implies that 
the wave Is subsonic and not a shock. However, 
we have applied an additional criterion, 
not discussed above, which Is not fulfilled by 
this otherwise acceptable solution. 


AV II n, AVl AB 



HOWEVER, RELATIONS ARE: 



Figure A. Relative orientations of the velocity 
and magnetic field Jumps. The upper 
panel shows the orientation antici- 
pated for a quasl-perpendlcular 
shock or wave, The bottom panel 
shows the actual orientations of the 
vectors used In the analysis presen- 
ted lr. Table 2 and which differ 
drastically from expectation. 


For both waves an d sh ocks, cer tain geometric re- 
lations between A V and A B must be satis- 
fied (Ref. 11, p. 98). In part icul ar, for a 
perpendicular shock or wave, A V s houl d be 
parallel to ft (and perpendicular to A B) as 
Indicated In the upper half of Figure A. The 
actual relation Is pres ente d In the lower half 
figure which shows that A V makes a large angle 
of ~ 72* with respect to n. As a consequ enc e , 
the numerators In V s , 1 .e ■ , A V • ft and A V x 
(?2, are small for a given A v. This mlsorlen- 
tatlon Is presumably the reason why V s Is low. 

The Ranklne-Hugonlot solutions (with 8= 10 and 
0 ■ 88*) lead to Jumps In density and field 
consistent with M ■ 1.6 + 0.3, l.e., a mean 

value similar to that obtained for the Inbound 
crossing. 


We are continuing to analyze the outbound cros- 
sing In the hope of Identifyin g t he origin of 
this mlsorlentatlon of A V and A B. It may be 
that time variations are responsible for pro- 
ducing this apparent discrepancy. The variabi- 
lity of the outbound crossing, which may be the 
result of multiple crossings of the bow wave or 
a consequence of the high value of B* 10, leads 
to a relatively long Interval of ~ 20 minutes 
between stable upstream and downstream condi- 
tions. It may be that the solar wind, specifi- 
cally the orientation of B, has changed signi- 
ficantly during this Interval. Continued analy- 
sis should show whether or not It Is possible 
to obtain an unambiguous result for the out- 
bound crossing. 
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